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Abstract 

Interferon-induced proteins with tetratricopeptide repeats (IFITs) are involved in tine protective response to viral infection, 
altliougli the precise mechanism of IFITs for reducing viral proliferation is currently unknown. The interaction with the 
translation initiation factor elF-3 or viral proteins and the sequestering of viral RNA have been proposed as potential 
antiviral functions for these proteins. In humans, four members of this family have been characterized. Nevertheless, 
information about these proteins in fish is almost non-existent. Exploiting the conservation of synteny between human and 
zebrafish genomes, we have identified ten members of the IFIT family located on four different chromosomes. The 
induction of these genes was examined both in vitro and in vivo after interferon (IFN) administration and rhabdovirus 
challenge. Whereas an induction of IFIT genes was observed after interferon treatments (IFNcDl, IFN<1)2 and IFN03), the viral 
infection did not affect these IFN-induced genes in vitro, and even reduced the IFN-induced expression of these genes. The 
response was largely different in vivo, with a broad up-regulation of IFIT genes after viral challenge. In addition, three 
selected IFITs were cloned in an expression vector and microinjected into zebrafish larvae to examine the protective effect 
of IFITs upon viral infection. Reduction in the mortality rate was observed confirming a conserved antiviral function in non- 
mammalian species. 
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Introduction 

Host antiviral innate immune responses begin witli the detection 
of viruses, which triggers the induction of cellular and molecular 
effectors with broad antiviral activities [1], including type I 
interferon (IFN) and hundreds of IFN-stimulated genes (ISGs), 
which contribute to the overall effects against a given virus [2,3] . 

In fish, although numerous IFN genes have been characterized 
(reviewed in [4]), their classification is controversial, as these genes 
are more diverse than previously thought, and their genomic 
structures, containing introns, are reminiscent of mammalian type 
III interferons (IFN-X), although the encoded proteins are similar 
to type I interferons [5]. Fish type I IFNs, named as interferon-phi 
(IFNfl>) were classified in two groups: group I (comprising IFNcDl 
and IFNm) and group II (composed by IFNm and IFN^S) [6,7] . 
These two groups of IFNs do not bind to the same receptor 
complexes in zebrafish, as was shown by Aggad et al. [8] . 

Knowledge of the antiviral properties of individual ISGs is 
mostly limited to a few intensively studied examples, such as PKR 
[9] or MX [10,11]. MX is a well-studied ISO in fish but our 
knowledge of other ISGs (e.g., Vig- 1 /viperin, ISG15, fiiiTRIMs, 



and PKR) in fish is limited (reviewed in [12]). Even in mammals, 
although the importance of the IFN system is clear, there is not 
enough information concerning the mechanism underlying the 
IFN-mediated inhibition of viral replication, particularly in vivo 
[13]. 

Among ISGs, a protein family called IFIT (interferon-induced 
proteins with tetratricopeptide repeats), which is characterized by 
tetratricopeptide repeats (TPR domains), has been examined in 
higher vertebrates [14—16]. Recent studies have shown that the 
IFIT family is conserved in mammals, amphibians and birds, but 
these proteins are not present in yeast, plants or lower animals, 
such as fruit fly and nematodes [16]. The members of this protein 
family were initially named according to their molecular weights 
(ISG54/P54, ISG56/P56, ISG58/P58 and ISG60/P60), although 
currently the most relevant and recent publications have adopted 
the term IFIT [16]. IFIT proteins are involved in many processes 
in response to viral infection and other functions, such as protein- 
protein and protein-RNA interactions, double-stranded RNA 
signaling, cell migration, and proliferation [14,17]. The transcrip- 
tional induction of the IFIT family genes has been described after 
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infection with both DNA- and RNA-viruses [18-21] and after 
bacterial stimulation in a type I IFN-dependent manner [22-24]. 
Although their antiviral mechanisms are still poorly understood, 
studies have shown that IFIT genes restrict virus replication 
through the alteration and suppression of protein synthesis or 
direct binding and sequestering of viral RNA, thereby reducing 
their infectivity [16]. 

In fish, information concerning IFIT genes is almost non- 
existent. Only a few partial or unconfirmed sequences of IFIT 
genes have been identified using sequencing analyses [16,25]. A 
formal characterization of IFIT genes and an in depth study of 
their regulation under different stimuli has never been done in fish. 

Teleosts fish offer an interesting model for the study of IFITs, 
not only for the clear interest in this ISG family in relation to viral 
infection, which constitutes an important threat, particularly for 
cultured fish, but also due to the ancient separation of fish from 
tetrapods and the great diversity of fish species. In addition, the 
advantage of the increased use of zebrafish [Danio rerio) as an 
important vertebrate model for studies in developmental and 
biomedical research, hematopoiesis and recently, immunology, 
has facilitated the development of genomic tools that allow the 
identification of new gene families. In the present work, we 
describe the complete repertoire of IFIT genes in zebrafish. Our 
study reveals a protein family forged through ancient duplication 
events, according with the results recently published [25]. To 
further explore the antiviral properties of these IFN-stimulated 
genes, in vivo and in vitro experiments were conducted in zebrafish 
after treatment with different recombinant IFNsfl> and after viral 
infection. Moreover, the protective effect of three selected 
zebrafish IFITs upon viral challenge was also examined in vivo. 

Results 

Defining the Complete Repertoire of IFIT Genes in 
Zebrafish 

Using a zebrafish genome-wide blast search, we detected a high 
degree of synteny between the human chromosome 10 (region 
q23.31) and the zebrafish chromosomes 17 and 12 (Figure lA). 
Our analysis confirmed the presence of five and three IFIT genes 
on zebrafish chromosomes 12 and 17, respectively. Moreover, 
another two genes were identified as similar to IFITs (one gene on 
chromosome 5 and the other gene on chromosome 13). 

The zebrafish IFIT sequences were named according to their 
chromosomal position: IFIT5A (ENSDARG00000088069), 
IFIT12A (ENSDARG00000008098), IFIT12B (EN- 
SDARG()()()()()0()7467), IFIT12C (ENSDARG00()00()90537), 
IFIT 1 2D (not identified in the Ensembl database), IFIT12E 
(ENSDARG00000090977), IFIT13A (ENSDARG00000057173), 
IFIT17A (ENSDARG00000071012), IFIT17B (EN- 
SDARG00000043584) and IFIT17C (ENSDARG000000,%976). 
In order to confirm the IFIT sequences, we designed specific 
primers to amplify and sequence the 10 zebrafish IFIT genes 
(primers in Table SI). The confirmed full-length ORFs were 
submitted to GenBank under accession numbers KF418356- 
KF418365. 

The block of human IFIT genes and the pseudogene IFIT IB 
(IFIT-1, IB, 2, 3 and ,5) is flanked downstream by the SLC16A12 
(solute carrier family 16, member 12) and PANKl (pantothenate 
kinase 1) genes. These genes are also present on zebrafish 
chromosomes 1 7 and 1 2 but are situated upstream of the IFIT 
region (Figure lA). The FAS (TNF receptor superfamily, member 
6) and CH25H (cholesterol 25-hydroxylase) genes are located 
upstream of the human IFIT cluster and showed sub-partitioning 
in zebrafish. Thus, FAS is located on chromosome 17 and CH25H 



is located on chromosome 12, and both genes have an inverted 
orientation at the 3' region of the IFIT genes (Figure lA). The 
q22.2 region of human chromosome 10 also showed homology 
with zebrafish chromosome 13. One IFIT-related gene (IFIT 13 A) 
was identified between the COMTDl (catechol-O-methyltrans- 
ferase domain containing 1) gene at the ,5' end and the NEFH 
(neurofilament, heavy polypeptide) and VDAC2 (voltage-depen- 
dent anion channel 2) genes at the 3' end of this chromosome. In 
this case, synteny was not conserved because there were no IFIT 
genes between the VDAC2 and COMTDl genes in human 
chromosome 10 (Figure lA). Moreover, another IFIT gciu- \vas 
identified on chromosome 5, but it was not possible to identify a 
conserved region between both species (Figure lA). 

Phylogenetic Tree and Analysis of Darwinian Selection 

Data mining for IFIT protein sequences retrieved 77 different 
sequences, 60 of which belonged to mammalian species, 2 bird 
species, 6 amphibian species and 9 fish species (Table S2). We used 
these sequences to construct a phylogenetic tree, revealing a clear 
separation of the IFIT sequences among the analyzed taxonomical 
classes (Figure IB). The tree topology categorized the mammalian 
IFIT sequences into four main branches, IFITl/lB, IFIT2, IFIT3 
and IFITS homology groups, with great confidence values; 
however, the phylogenetic relationship among non-mammalian 
IFIT sequences was not clear, potentially reflecting the minor 
representation of these sequences in the entire analysis. For 
example, the evolutionary relationship of the different IFIT genes 
across fish species was not confidently resolved, although some 
sequences were orthologous. However, an interesting pattern 
emerged, branching the zebrafish IFIT sequences belonging to the 
same chromosome. 

We estimated the dN/dS ratios ((o) among the zebrafish IFIT 
sequences to quantify the selection pressure acting on IFITs genes 
and determined that the genes located on chromosome 1 7 
underwent positive Darwinian selection (co>l). The dN/dS ratio 
observed between IFIT 1 7 A and IFIT 1 7B was 1 .2928, whereas this 
value was higher between IFITs 17B and 17C ((0 = 1.4667) and 
between IFITs 17A and 17C (m = 1.5130). The IFIT genes located 
on chromosomes 5, 12 and 13 were not subjected to this 
evolutionary mechanism, obtaining CO values lower than 1 . 

Sequencing Analysis and Structure Domains 

The study of the domain structure revealed the presence of the 
characteristic TPR motifs in all the analyzed sequences (Figure 2A), 
but the number of these repeats was variable, ranging from two up 
to eleven repetitions. However, we also observed variability in the 
number and position of these characteristic domains in the four 
human genes and in the three IFITs described in mice (Figure 2A). 
The values of identity and similarity between zebrafish and human 
and murine proteins were lower than 37% and 62%, respectively 
(Table S3). The number of amino acids varied between 302 and 
483 residues for all sequences, except for isoform 5 A, which 
presented a total of 1038 amino acids. Thus, differences in 
calculated molectilar weights were also observed, whereas most of 
the proteins presented molecular weights ranging from 50- 
56 kDa, tiie 12D isoform was 34.92 kDa and IFIT5A was 
120.29 kDa (Figure 2B). Regarding the theoretical isoelectric 
point (pi), the isoforms located on chromosomes 5, 12 and 17 
presented values low(;r than 7.0, except 12A, whose pi value was 
8.20. IFIT13A also showed a more basic pi value of 8.75 
(Figure. 2B). 

We also examined the tridimensional structure of the IFITs and 
identified three different structural models (Figure 20). The TM- 
scores observed for zebrafish IFITs revealed that human IFIT5 
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Figure 1. Synteny and phylogenetic analysis. A. Comparative gene location of IFIT sequences betw/een human chromosome 10 and zebrafish 
chromosomes 5, 12, 13 and 17. The genes flanking the IFIT cluster, located at the q23.31 region of human chromosome 10, shovj/ed a high degree of 
synteny with the zebrafish chromosomes 17 and 12. Moreover, the genes surrounding zebrafish IFIT 13 were located on the same human 
chromosome but in the q22.2 region. An IFIT sequence was also present on zebrafish chromosome 5. B. Phylogenetic relationship between the IFIT 
proteins from bony fish and similar proteins described in other vertebrates, such as amphibians, birds and mammals. Lch: Latimeria chalumnae; Om: 
Oncorhynchus mykiss; Sas: Salmo salar; Ca: Carassius auratus; Ci: Ctenopharyngodon idella; Lc: Larimichthys crocea; Dl: Dicentrarchus labrax; Tn: 
Tetraodon nigroviridis. 
doi:1 0.1 371 /journal.pone.01 0001 S.gOOl 

was the main analog protein for IFIT5A (TM-score = 0.448), the superhelical TPR-repeat domain of O-linked GlcNAc 
IFIT12A (TM-score = 0.981), IFIT12B (TM-score = 0.966), IFIT transferase was the template for the construction of the IFIT 
12C (TM-score = 0.951), IFIT 1 2D (TM-score = 0.950), IFIT12E 17C SD-structure (TM-score = 0.813). 
(TM-score = 0.969) and IFIT13A (TM-score = 0.952). However, 
human ISG54 or IFIT2 was the best analog for IFIT17A (TM- 
score = 0.947) and IFIT17B (TM-score = 0.943), and interestingly. 
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Figure 2. Sequence analysis and structural domain of the zebrafish IFITs. A. Tetratricopeptide repeat (TPR) motifs on the 10 IFIT proteins 
present in zebrafish. Orthologous proteins from human and mouse are shown for comparative purposes. B. Length, calculated molecular weight 
(kDa) and theoretical pl of the different zebrafish IFITs. C. 3D-structure of zebrafish IFITs, predicted using 1-TASSER server, selecting the model with the 
best C-score and viewed using PyMOL 
doi:1 0.1 371 /journal.pone.01 0001 5.g002 



Constitutive and Tissue-specific Expression of IFIT Genes 

The analysis of 8 different adult zebrafish tissues revealed a 
higher basal expression of IFIT genes on chromosome 12, being 
IFIT12C the gene with the largest presence in the whole of the 
tissues analyzed (spleen, kidney, gill, caudal fin and head). 
Interestingly, IFITs with low expression in most of the tissues, 
showed significantly higher expression in the intestine (12A, 17 A, 
17B) or muscle (13A, 17C) (Figure 3A). Regarding to the relative 
proportion of the IFIT genes in the analyzed tissues, whereas in 
spleen, kidney, muscle, intestine and liver all the IFITs were 
present, in giUs IFITs from chromosome 12 were predominant 
(Figure 3B). 

The constitutive expression of IFIT genes was also analyzed in 
ZF4 and kidney primary cells (Figure SI). ZF4 cells presented 
higher expression levels than kidney primary cell cultures. In 
addition, IFITs from zebrafish chromosome 12 showed a higher 
basal expression in both cell types than those from chromosomes 



5, 13 or 17, except for isoform 12B, which showed lower 
expression in ZF4 cells. 

Interferons Induce the Expression of IFIT Genes Inlvitro 

First, we determined whether three selected zebrafish interfer- 
ons (IFN<1)1, IFN02 and IFN<1)3) induced the expression of 
zebrafish IFITs. Therefore, we analyzed the biological activity of 
recombinant zebrafish IFNsCD. The activity of supernatant from 
HEK-293 cells transfected with plasmids containing IFNs 
sequences was first confirmed by a decrease in the viral titer of 
spring viraemia of carp virus (SVC V) and the induction of MXab 
(isoforms a and b) expression in ZF4 and kidney primary cell 
culture, as shown in the Figure S2A and S2B. 

The treatment of ZF4 cells and kidney cell cultures with 
interferons induced changes in the expression of IFIT genes. 
Overall, the results showed higher expression values for all IFIT 
genes in kidney cells than in ZF4 cells (Figure 4A and 4B). As 
expected, these results suggest a role of IFNs in the induction of 
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Figure 3. Tissue-specific expression of zebrafish IFIT genes. A. Constitutive expression of IFIT genes in tissues of adult zebrafish (S: Spleen; L: 
Liver; K: Kidney; G: Gill; CF: Caudal fin; H: Head; M: Muscle; I: Intestine). For basal expression of each IFIT form, tissues were sampled and pooled, 
yielding a total of 4 pools of 5 fish per organ. The relative expression level of each gene, normalized to the expression level of the 1 8 S ribosomal RNA 
gene in the same tissue, was expressed as arbitrary units. The graphs represent the mean ± standard error of 4 independent samples. B. Relative 
proportion of the IFIT transcripts in different zebrafish tissues. 
doi:1 0.1 371 /journal.pone.01 0001 5.g003 
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Figure 4. In vitro effect of IFNs on IFITs expression in ZF4 cells and Itidney primary cell cultures. A. Expression levels of IFIT genes in ZF4 
cells after 4 and 24 hours of stimulation with supernatants from transfected HEK-293 cells with plasmids containing sequences for zf-IFN<I)1, zf-IFNil>2 
and zf-IFNOS. B. Kidney primary cell cultures after 4 and 24 hours of stimulation with supernatants from transfected HEK-293 cells with plasmids 
containing sequences for zf-IFNOI, zf-IFN<I)2 and zf-IFNOB. After 4 and 24 hours of stimulation, RNA was extracted and the cDNA was synthesized. 
Analysis of the gene expression was performed through real-time PCR, using 18 S ribosomal RNA as a housekeeping gene. The expression level of 
each gene was expressed as fold-change with respect to the empty plasmid group. The data are represented as the mean ± standard error of three 
independent samples. Significant differences among cells transfected with IFN<I) plasmids and empty plasmid were displayed as ***(0.0001<p< 
0.001), **(0.001<p<0.01) or *(0.01<p<0.05). 
doi:1 0.1 371 /journal.pone.01 0001 5.g004 



zebrafLsh IFITs. Moreover, the fact that IFITs basal expression 
was higher in ZF4 than in primary cell cultures and that these cells 
showed a higher response to IFNs suggest a cell-specific response 
or mechanism. 

In ZF4 cells, treatment with IFNs after 4 hours did not induce 
high levels of IFITs expression, with the exception of IFIT17C, 
which was significandy induced, regardless of the IFNfl) used. 
Interestingly, a significant down-regulation of isoforms 5A, 1 2A, 
12C, 12D and 13A was observed in cells treated with the IFNOS 
(Figure 4A). A higher number of isoforms was significandy induced 
at 24 hours compared with the results observed at 4 hours, and the 
IFIT on chromosome 5 was the most induced through IFN(I>2 and 
IFN<I)3. The IFITs on chromosome 12 (12B, 12C and 12E) were 
significantly induced, regardless of the IFN® used (Figure 4A). 

In primary cell culture, the IFIT genes located on chromosomes 
1 3 and 1 7 were more induced than those on chromosomes 5 and 
12, and isoform 13A showed the strongest induction at 4 hours 
(Figure 4B). Regarding the effect of the different IFNs®, most of 
the IFITs on chromosomes 12 and 13 were induced through 
IFN$2 and IFN<I)3 at 4 hours (Figure 4B). Thus, in primary cell 
culture, only IFIT12B, IFIT12C, IFIT13A, IFIT17A and 
IFIT17C were significantly induced through IFNOl after 4 hours, 
whereas, after 24 hours, the significant induction of 12B, 12C, 
12E, 13A, 17A, 17B and 17C was observed. 



IVlodulation of IFITs Expression upon Viral Challenge 

Once we determined that zebrafish IFITs were modulated 
through interferons and showed different tissue expression profiles, 
we evaluated the effect of an in vitro viral infection on the 
expression of IFITs in ZF4 and kidney primary cells. Surprisingly, 
SVCV did not modify the expression of the different IFITs (with 
the exception of IFIT,5A in ZF4 and IFIT12B in kidney cells that 
showed a slight expression increase) (Figure 5 A). To determine 
whether this effect was induced through a direct effect of the virus 
on IFITs expression or if the virus was affecting the interferon 
signaling cascade, we measured the expression levels of IFNOl, 2 
and 3 and the interferon-induced protein MXab after in vitro 
infection. The results showed that IFN^tl, 2 and 3 and MXab 
were not induced in either ZF4 or kidney cells through a 24 hours 
in vitro infection (Figure 5B). 

Next, we determined whether the IFNO-induced expression of 
IFIT genes was also modulated through viral infection. Thus, the 
effect of the virus on kidney primary cells treated with the 
recombinant IFNs® was analyzed at 24 hours after infection 
(Figwe 6). In general, the virus reduced the expression of 
interferon-induced IFITs. However, the expression induced 
through recombinant IFNOl did not show this clear decrease, 
and in the case of isoforms 1 2A, an increment in expression was 
observed after viral infection. In addition, isoform 12A was the 
only gene that experienced an up-modulation after stimulation 
with the three recombinant IFNsO and the infection with the virus 
(Figure 6). 
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Figure 5. In vitro effect of viral infection in ZF4 cells and kidney primary cell cultures. A. Expression of IFIT genes in ZF4 cells and kidney 
primary cell culture at 24 hours after infection with SVCV. B. MXab and IFN<I>1, 2 and 3 expression in ZF4 cells and kidney primary cell culture at 24 
hours after infection with SVCV. After 24 hours of stimulation, the RNA was extracted, and the cDNA was synthesized. The analysis of gene expression 
was performed through real-time PCR, using 18 S ribosomal RNA as housekeeping gene. The expression level of each gene was expressed as fold- 
change with respect to the control group, non-infected cells. The data are shown as the mean ± error of three independent biological samples. 
Significant differences among infected and uninfected group were displayed as ***(0.0001<p<0.001), **(0.001<p<0.01) or *(0.01<p<0.05). 
doi:1 0.1 371/journal.pone.01 0001 5.g005 



To determine whether the trend of down-modulation of IFN- 
induced gene.s through SVCV infertion also occurred in vivo, we 
examined IFIT genes expression in kidney cells from adult animals 
injected with the virus. In this case, all the IFITs showed an 
increase in expression after 24 hours (the increased expression of 
5A and 17B were not statistically significant), and 17A showed the 
highest fold change (Figure 7A). As expected, an up-regulation of 
all analyzed IFNcDs was observed, being the expression of IFNCD2 
the highest detected. The interferon inducible protein MXab 
showed a statistically significant increase of 40 fold after the 
stimulation with SVCV (Figure 7B). These results denoted the 
different response to the virus after an inPvitro or inPvivo infection. 

Evaluation of Antiviral Activity of IFIT Genes in Zebrafish 
Larvae 

Next, we analyzed the in vivo antiviral activity of selected IFITs 
in zebrafish larvae previously microinjected with expression 
vectors containing either IFIT12B, 13A or 17 A, f()llo\\Td by 
infection with SVCV. When zebrafish larvae were infected with 
SVCV (control group), the mortalities reached the maximum level 
at 36 hours after challenge, and only 10% of the animals survived 
the infection (Figure 8A). At this point, the percentages of survival 
in animals treated with expression vectors containing IFIT 
sequences were higher than those obtained in the control group. 
Animals treated with 12B, ISA and 17A showed a final % survival 



of 17.5, 22.5 and 40%, respectively. At 36?h post-infection, 
significant differences in the survival were observed in animals 
treated with the three IFITs, whereas at 48 and 72 hours only the 
larvae treated with IFIT17A showed a significant increase in 
survival with respect to the control group. 

The transcription of the SVCV N gene was also measured 
through qPCR at 9 h after infection to determine the effect of 
these IFITs on viral transcription. At this time point, only larvae 
injected with the IFIT 1 7 A plasmid showed a significant lower viral 
N gene transcription compared with the infected control group 
(Figure 8B). 

Discussion 

IFlTs are a ncnel IFN-stimulated gene family with antiviral 
properties not formally described in fish until a recent published 
work [25]. Using genome synteny and sequence comparison, we 
identified 10 sequences in the zebrafish genome with homology to 
human IFIT genes located in chromosome 10 (four genes and one 
pseudogene). With the exception of IF1T5A, all zebrafish IFIT 
genes presented similar length and conserved domains. 

In contrast to humans, most of the IFIT family genes in 
zebrafish are located in two chromosomes (12 and 17), similar to 
the structure observed in dogs [14,16]. Five IFIT genes were 
clustered on chromosome 12, as previously described [16,25], but 
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Figure 6. /n i7fro effect of viral infection and IFN treatment in icidney primary cell cultures. Expression level of IFIT genes in kidney primary 
cell cultures after 24 hours of stimulation with supernatants from transfected HEK-293 cells with plasmids containing sequences for zf-IFNflJl, zf- 
IFN<I)2 and zf-IFN<I)3 in combination with SVCV. After 24 hours of stimulation, the RNA was extracted, and the cDNA synthesized. The analysis of gene 
expression was performed through real-time PCR, using 18 S ribosomal RNA as a housekeeping gene. The effect of the virus infection on the 
expression of IFITs induced by the different IFNs was represented as a fold-change with respect to the group stimulated with supernatant from cells 
transfected with the empty plasmid. The data are represented as the mean ± standard error of three independent samples. The asterisks denote 
significant differences between infected and non-infected groups. Significant differences were displayed as ***(0.0001 <p<0.001 ), **(0.001 <p<0.01 ) 
or *(0.01<p<0.05). 

doi:1 0.1 371 /journal.pone.01 0001 5.g006 
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Figure 7. In vivo effect of viral infection on IFNs, MXab and IFIT genes expression in Itidney. A. Expression of IFIT genes in kidney cells 
from adult zebrafish at 24 hours after infection with SVCV. B. Expression of IFNs and MXab in kidney cells from adult zebrafish at 24 hours after 
infection with SVCV. Adult individuals were injected intraperitoneally with 10 |J.I of SVCV (2.7x10^ TCIDjo/ml). RNA was isolated from head kidney 
cells, 24 hours post-infection. cDNA was obtained, and real-time PCR was performed using 18 S ribosomal RNA as a housekeeping gene. The 
expression level of each gene was expressed as fold-change with respect to the levels detected in the control group (injected with culture medium). 
The data are shown as the mean ± standard error of three individuals. The asterisks denote statistically significant differences with respect to the 
control group. Significant differences were displayed as ***(0.0001<p<0.001), **(0.001<p<0.01) or *(0.01<p<0.05). 
doi:1 0.1 371 /journal.pone.01 0001 5.g007 



our analysis also registered the presence of three IFIT genes on 
chromosome 17, one gene on chromosome 13 and another 
additional gene on chromosome 5, in agreement with that recently 



reported [25]. The IFIT genes on chromosomes 5 and 13 did not 
conserve the synteny with vertebrate chromosomes, likely reflect- 
ing genomic translocations [26,27]. Genome duplication events 
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Figure 8. Antiviral activity of IFITs 12B, 13A and 17A in larvae infected with SVCV. A. Antiviral activity of selected IFITs was evaluated in 
zebrafish larvae. One-cell stage zebrafish embryos were microinjected with 100 pg/egg (final volume of 2 nl) of the recombinant plasmids pcDNA 
3.1-IFIT12B, pcDNA 3.1-IFIT13A, pcDNA 3.1-IFIT17A as well as pcDNA 3.1 -empty. Three days after plasmid administration, the larvae were 
microinjected in the duct of Cuvier with 2 nL of a SVCV suspension at a final concentration of 10^ TCIDso/ml. The data are shown as the percentage of 
survival observed at 3 days after infection. Significant differences (PsO.05) in the percentage of survival between larvae treated with the IFITs and the 
control group are indicated with asterisks. The results are represented as the mean ± standard error of four independent samples. B. The relative 
expression level of the viral N gene was analyzed through qPCR at 9 hours post-challenge. The raw data were normalized using the 18 S ribosomal 
RNA as a housekeeping gene. The results are presented as the mean ± standard error of three biological replicates. Significant differences were 
displayed as ***(0.0001<p<0.001), **(0.001<p<0.01) or *(0.01<p<0.05). 
doi:1 0.1 371/Journal.pone.01 0001 5.g008 
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are powerful drivers of evolution [28,29], as they provide 
opportunities for the modification or mutation of the gene 
duplicates, while critical functions are maintained through the 
other copies. The modified or mutated duplicated genes might 
acquire new functions or divide the original functions of the 
ancestral gene among diflFerent isoforms. The presence of difiFerent 
members of the IFIT family genes in fish might facilitate the 
expansion of innate immune recognition or modulate innate and 
adaptive immune responses to specific challenges. 

The analysis of Darwinian scdection was conducted to quantify 
the selection pressures acting on IFIT genes [30,31], and the 
results showed that the genes located on chromosome 17 
underwent positive Darwinian selection. This result, together with 
the phylogenetic analysis, reflected the increased accumulation of 
evolutionary changes in these genes with respect to the other IFIT 
genes. The accelerated evolution of the IFIT genes on chromo- 
some 1 7 might be associated with the direct interaction of these 
proteins with pathogenic viruses, as an elevated selective pressure 
and rapid evolution of immune-related genes, particularly those 
that directly interact with pathogens, was observed compared with 
non-immune genes [32-41]. 

Phylogenetic analysis confirmed that mammal IFITs are 
clustered in a main group comprising the four IFIT g(;n(;s 
previously described [14,16]. The sequences from other verte- 
brates constituted diflFerent clusters, depending on the class 
(amphibians, birds or bony fish), as previously reported [16] 
however, the construction of an internal classification of the IFIT 
genes in fish is difficult due to the scarce information available in 
public databases for other fish species. A deeper analysis of this 
gene family for organisms belonging to different taxa would help 
elucidate the evolutionary process involved. 

The functional activity of the novel IFIT family described in 
zebrafish was explored using in vitro and in vivo experimental 
models. Most mammalian cell types do not express IFIT genes 
under basal conditions [15,42]; however, \vc observed a constitu- 
tive expression of the IFITs particularly the genes located on 
chromosome 12. In addition, the distribution of the ten diflFerent 
IFIT transcripts in zebrafish displayed distinct patterns of 
preferential expression at tissue level, and revealed an extremely 
higher functional complexity than that previously reported in mice 
[14,20,43,44]. Constitutively expressed genes, such as 12C, 12D 
and 12E, were present in all tissues, but the low expression of other 
genes, such as 1 2B, 1 3A and 1 7A, might suggest they are inducible 
genes. Moreover, the high expression levels observed in the 
intestine might suggest a specific function for IFIT12A, 17A and 
17B in this tissue. Indeed, the IFITs on chromosome 17 are 
primarily expressed in the liver and intestine. These results suggest 
that IFIT genes might have non-redundant antiviral functions, as 
previously suggested in mice [20,44], reflecting the differentiation 
and subsequent specialization of the members, potentially 
facilitated through the gene expansion observed in fish. 

The expression of IFIT genes was analyzed in kidney primary 
cell cultures and in ZF4 cells in response to an IFNsO treatment. 
Recombinant zebrafish IFNs€> (1, 2 and 3) significantiy reduced 
the viral titer in ZF4 cells infected with SVCV, as previously 
reported [4,5], and induced a rapid and high expression of the 
interferon-inducible MXab genes in kidney primary cell cultures 
and ZF4 cells, as described in other fish models. Interestingly, the 
IFNsfl> from group II (IFN<I)2 and 3) showed higher antiviral 
activity and induction of the MXab genes than that observed for 
IFNsO from group I (IFN<1>1). This differential antiviral activity 
observed between IFNs from group I (IFNcDl) and group II 
(IFN<I)2 and 3) could reflect the induction of several response 
pathways, as these molecules do not bind the same receptor 



complexes [8]. The treatment of cells with IFNsO also induced a 
rapid increase in IFIT genes expression in kidney cells (mainly 
IFIT13A and 17s) consistent with previous studies [16]. The 
modulation of IFIT genes in ZF4 cells was much lower than that 
observed in kidney cells most likely because ZF4 is not an immune 
cell line [46], and the effect of IFN stimulation was not 
comparable with the effect observed in specific immune cells 
presented in kidney primary cell cultures and the hematopoietic 
tissues of the fish. 

During evolution, some viruses have evolved sophisticated 
mechanisms to avoid the host innate immune system. In 
particular, some rhabdoviruses, such as human virus VSV 
(vesicular stomatitis virus) or RV (rabies virus), have developed 
counteractions to both IFN induction and IFN signaling [47]. 
These viruses have diflFerent mechanisms for antagonizing the type 
I interferon response and blocking the induction of antiviral 
molecules; however, in both cases, the objective is the evasion of 
the host immune defense [48] . In fish, the matrix protein of the 
novirhabdovirus, IHNV (Infectious hematopoietic necrosis virus), 
affects host celMar gene expression to inhibit the transcription of 
immune-related genes [49]; however, littie is known about how 
this effect is orchestrated. 

ZF4 cells and kidney primary cell ctdtures infected with SVCV 
did not show a typical anti-viral response upon IFN gene induction 
and the increased expression of ISGs, such as MX or IFITs [2,50]. 
The blocking of the interferon system suggests that the virus 
suppresses the immune response in primary celt cultures.This 
response was also investigated when kidney cells were forced to 
mount an antiviral response through the stimulation with 
recombinant IFNsO and also were infected with the virus SVCV. 
In this case, the IFITs expression pattern was modulated, as 
described in human hepatocytes infected with hepatitis C virus 
[51]. Kidney cells treated with IFN<1)2 and 3 showed the reduced 
expression of almost all IFIT genes (IFIT 1 2 A made the difference), 
whereas cells treated with IFNfDl only showed the down 
modulation of IFIT17C after viral infection. This result could 
indicate that the virus avoids the host defense system and 
suppresses the expression of a specific subset of IFIT genes for 
the establishment of infection. The viral-mediated inhibition of the 
IFN system has been previously described [51,52,53]. The 
response pattern observed in cells treated with IFNcDl after 
SVCV infection might reflect different signaling pathways between 
cells stimulated with IFNsO from groups I and II. Importantly, 
IFIT12A was the only gene whose expression was synergistically 
induced through all interferons and in response to virus infection. 
The different behavior after in vitro stimulation together with the 
tissue-specific expression of IFITs genes suggests the expansion 
and differential functions of these genes. 

However, when viral infection is conducted in vivo, after 
intraperitoneal infection, a clear up-regulation of the expression 
of IFNs and ISGs, including MXab and IFITs, was observed. The 
overall immunity of the host is required to orchestrate the effective 
control of viral infections, and the absence of a complete response 
often results in fatal infections [52]. Both ZF4 cells and primary 
ceU head kidney leukocyte cultures exhibited limited defense 
against viral infections because of the incomplete host machinery 
and these models are, therefore, easily manipulated by SVCV. 
These cells respond to IFN stimulus but are unable to mount an 
effective response against virus. 

The most induced IFITs after in vivo infection, 12B, ISA and 
17A, were selected to confirm direct antiviral activity. The 
microinjection of zebrafish eggs at one-cell stage with expression 
plasmids encoding these genes induced a significant reduction in 
mortality after SVCV infection, highlighting the antiviral role that 
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these proteins might play in non-mammaUan species. In mam- 
mals, there is evidence implicating these proteins in the restriction 
of translation initiation through interactions with the translation 

initiation factor eIF-3 [44,54-57]. Moreover, IFITs are able to 
sequester viral proteins, such as human papillomavirus helk:ase El 
[58] and inhibit virus replication through the direct binding and 
sequestering of viral nucleic acids [59-63]. However, it remains 
unknown whether the same mechanisms are also present in fish. 
We can confirm that the present results provide the basis for 
multiple future research studies concerning not only the protection 
of fish (particularly aquacultured species) against virus infections 
but also the investigation of the basic aspects of IFITs biology, 
which could be studied in zebrafish, an attractive model organism 
with numerous experimental advantages. 

Materials and Methods 

Sequence Retrieval and Analysis 

The IFIT sequences were searched using the zebrafish genome 
assembly version Zv9 (www.ensembl.org/Daniorerio/), exploiting 
the conservation of synteny between the human and zebrafish 
genomes. The sequences were confirmed through PGR amplifi- 
cation using specific primers (Table SI) to obtain the full-length 
open reading frame (ORF) of each gene. The PGR products were 
subcloned into a pCRS. 1 vector (Invitrogen) and transformed into 
One Shot TOPI OF' competent cells (Invitrogen) for subsequent 
sequencing and ORF confirmation. 

The identity and similarity analysis between the zebrafish, 
human and mouse IFIT sequences was performed using MatGAT 
[64]. The TPR distribution was analyzed using TPRpred (http:// 
tprpred.tuebingen.mpg.de/) [65], and the theoretical isoelectric 
point (pi) and the calculated molecular weight were determined 
using ExPaSy tools (http:/ /us. expasy.org/ tools). The 3D-structure 
of zebrafish IFITs was predicted using the I-TASSER server [66], 
selecting the model with the best C-score, and viewed through 
PyMOL (http://www.pymol.org). The Template Modeling Score 
(TM-score), a measure of structural similarity between two 
proteins, was also analyzed to identify structural analogs with 
known crystal architecture in the Protein Data Bank (PDB; http:// 
www.rcsb.org/ pdb/). 

Phylogenetic Tree and Analysis of Darwinian Selection 

IFIT-famUy protein sequences were r(;trie\cd from the NGBI 
Protein, Uniprot and Ensembl databases based on annotation. 
The sequences were subsequently complemented using a blastp 
search for homologs in different databases. The initial sequence 
alignment was performed using the MAFFT online server 
following an E-INS-i strategy [67]. The resulting alignment was 
pruned using Gblocks 0.91b [68] and subsequentiy analyzed using 
ProtTest 3.2 [69] to determine the best-fit amino acid replacement 
model using the Akaike Information Griterion (AIG) [70], specified 
to estimate the maximum likelihood gene tree using PhyML 3.0 
[71]. The nodal confidence \\-as calculated using the aLRT 
method [72]. Edition and representation of the obtained tree was 
performed in FigTree vl.3.1 (http://tree.bio.ed.ac.uk/software/ 
figtree/). 

An estimation of the rates of synonymous (silent) and 
nonsynonymous (amino-acid-changing) substitutions was p(;r- 
formed to identify a positive Darwinian selection in the zebrafish 
IFIT family using the PAML package version 4 [73]. The 
maximum-likelihood (ML) approach was implemented in the 
CODEML program to determine the CO value among zebrafish 
IFITs. 



Animals 

Wild type adult zebrafish [Danio rerio) were grown in our 
experim("ntal facilities according to established protocols [74,75] 
(also see iittp://zfm.org/zf_info/zfbook/ zfbk.html). Fish care and 
the challenge experiments were conducted according the GSIG 
National Committee on Bioethics under approval number 
07_09032012. 

Cell Cultures and Viral Infection 

Fibroblastic like cell line, ZF4, derived from 1 -day-old zebrafish 
embryos (ATCC CRL-2050) [46] were cultured in Dulbecco's 
modified Eagle's medium (D/MEM/F12, Gibco) supplemented 
with 100 |J,g/mL of primocin (InvivoGen) and 10% fetal bovine 
serum (FBS) at 26°C. Human HEK-293 cells (ATGG GRL-1573) 
[76] were grown in Eagle's Minimum Essential Medium (Gibco) 
supplemented with 100 |i.g/mL primocin (InvivoGen), IX non- 
essential amino acids (Gibco), 1 mM sodium pyruvate (Gibco) and 
10% FBS. The cells were incubated in a 5% CO2 atmosphere at 
37°C. Kidney cell suspensions were obtained from adult fish 
sacrificed using anaesthesia in ice. Kidneys were homogenized 
through a lOO-fim mesh, and the mixture was adjusted to the 
required concentration (1.5x10^ cells/ml) in Leibovitz L-15 
medium (Gibco) supplemented with 100 Hg/mL of Primocin 
(InvivoGen) and 2% FBS and maintained at 26°G. For the in vitro 
stimulations, the cells were seeded into 24-well plates at 1 ml per 
well. 

The rhabdovirus, spring viraemia of carp virus (SVGV isolate 
56/70) was used in these experiments. Experimental infections 
were performed at 22°G, and the viral titer was calculated as 
previously described [77]. 

RNA Extraction and Gene Expression 

Total RNA isolation was performed using the Maxwell 16 LEV 
Simply RNA Tissue Kit (Promega, Madison, WI, USA) according 
to the manufacturer's instructions. cDNA was obtained from 1 |J,g 
of total RNA using the SuperScript III First-Strand Synthesis 
Supermix (Invitrogen). Specific qPGR primers were designed 
(Table 81) using the PrimerS program [78], and the primer 
efficiency was evaluated [79]. A previously described [80] cDNA 
template was used for real-time PGR amplification, with 40 cycles 
and a 60°G annealing temperature. All reactions were performed 
with several biological replicates and using technical triplicates. 
The relative expression levels of the genes were normalized to the 
expression of 18 S ribosomal RNA [81], as a housekeeping gene 
control (primers specified in Table 81), following the PfafH method 
[79]. 

Production of Zebrafish Recombinant IFNs(p 

Zebrafish imm, IFN02 and IFN03 (GenBank accession 
numbers: NM_207640, NC_0071 14 and NC_0071 14, respective- 
ly) expression constructs in the pcDNA3.1/V5-His backbone were 
kindly provided by Dr. Mulero (University of Murcia, Spain). 
Recombinant IFNsfl) were produced by transfection of 6 |Xg the 
plasmids into HEK-293 cells at 70-80% confluence using the X- 
tremeGENE HP DNA Transfection Reagent (Roche) according to 
the manufacturer's instructions. Forty-eight hours after transfec- 
tion, the supernatants were collected and stored at — 80°G until 
further use. 

Antiviral Activity of IFIT Genes in Zebrafish 

Three selected IFITs (12B, 13A and 17A) were amplified using 
touchdown PGR (primers in Table SI), and the PGR products 
were cloned using the pcDNA 3.1/V5-His TOPO TA Expression 
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Kit (Invitrogen). One Shot TOPI OF competent cells (Invitrogen) 
were transformed to generate the plasmid constructs. Plasmid 
purifications were conducted using the Purelink HiPure Plasmid 
Midiprep Kit (Invitrogen). The recombinant pla.smids were 
microinjected into one-cell stage zebrafish embryos with a glass 
microneedle using Narishige MN-15I micromanipulator and 
Narishige IM-30 microinjector. In each experiment, a total of 
240 embryos were divided into 6 groups of 40 eggs (4 replicates of 
10 embryos) and each batch was microinjected with the following 
treatments dfluted in PBS: pcDNA 3.1 -IFIT 12B, pcDNA 3.1- 
IFIT13A, pcDNA 3.1-IFIT17A, pcDNA 3.1-empty, and PBS. An 
additional untreated group was included to control the egg quality 
and survival. The quantity of plasmid inoculated into each embryo 
was 100 pg/egg in a final volume of 2 nL. Three days after 
plasmid administration, the larvae were microinjected in the duct 
of Cuvier with 2 nL of a S VC V suspension at a final concentration 
of 1 0'^ TCIDjo/nil. The mortalities due to the viral infection were 
registered for 3 days after infection before an independendy 
feeding and therefore before an ethical approval is required (EU 
directive 2010_63) [82]. Fish condition was controlled three times 
a day. The viral transcription in IFIT-injected larvae was 
quantified through qPCR using specific primers for the N gene 
of the SVCV at 9 hours after the infection. The relative expression 
of the N gene was normalized to the expression of 18 S ribosomal 
RNA (primers specified in Table SI). 

Statistical Analysis 

The results were expressed as the means ± standard error. The 
significant differences were determined using Students t-test. The 
data from the in vivo antiviral activity of IFITs was analyzed using 
one-way analysis of variance (ANOVA) followed by Tukeys 
multiple comparison test. 

Supporting Information 

Figure SI Constitutive expression of IFIT genes in ZF4 
cells and in head kidney primary cell cultures. The basal 
expression of the different IFIT genes was analyzed through real- 
time PGR in ZF4 cells and in leukocyte primary cell cultures from 
kidney. The relative expression level of the genes was normalized 
using the 18 S ribosomal RNA as a housekeeping gene. The 
graphs represent the mean ± standard error of three independent 
samples. 
(TIF) 

Figure S2 Biological activity of recombinant zebrafish 

IFNs. A. The biological activity of the supernatants from HEK- 
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